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ABSTRACT 



Since the cosmic peculiar velocity field depends on small wave-number modes strongly, we 
cannot probe its universal properties unless we observe a sufficiently large region. We calculate the 
expected deviation (sample variance) of the peculiar velocity dispersion from its universal value in 
<^ , the case observed volume is finite. Using linear theory we show that the sample variance remains 
as large as ~ 10%, even if the observed region is as deep as 100/i _1 Mpc and that it seriously affects 
OQ | the estimation of cosmological parameters from the peculiar velocity field. 

Subject headings: cosmology: theory — large-scale structure of the Universe — methods: statistical 

^ ■ 

Q ■ 1. Introduction 

& : 

Observational analysis of the large-scale peculiar velocity field is considered as a very effective 
method to impose constraint on the cosmological parameters as well as the spectrum of the 
primordial density fluctuations in the universe. Since its power spectrum depends more strongly on 
small wave- number modes, which are less contaminated by nonlinear effects, than that of density 
fluctuations, linear perturbation theory suffices to reproduce the amplitude of the peculiar velocity 
even on a relatively small length scales. In fact, Bahcall et al. (1994) have calculated the one-point 
peculiar-velocity dispersion smoothed with a Gaussian filter over the smoothing scale 3/i _1 Mpc 
in their large-scale N— body simulations with various cosmological models and have found an 
excellent agreement between their calculations and the predictions of linear theory. Here h is as 
usual the present Hubble parameter in unit of lOOkm/sec/Mpc. Validity of the linear analysis has 
also been confirmed analytically in the framework of higher order Eulerian perturbation theory 
by Makino, Sasaki & Suto (1992), who showed that second-order effects were negligible on the 
smoothed velocity dispersion in the case of cold-dark-matter power spectrum. 

The stronger dependence on the smaller wave-number modes, on the other hand, implies that 
the correlation length of peculiar velocity field is larger than that of the density field and that we 
must survey a volume large enough to contain sufficient number of such modes. In that sense, 
Bahcall et al. (1994) was quite correct in adopting a simulation box whose dimension was as large 
as 800/i _1 Mpc. At present, however, we can by no means hope to observe the peculiar velocity 
field that far in the real Universe. For example, the observed depth of the recent high quality 
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catalogue by Giovanelli et al. (1997) has been limited to 100fo _1 Mpc mainly due to the difficulties 
to estimate the distance which is of course essential in obtaining the line-of-sight peculiar velocity 
(see also Giovanelli et al. 1996). Thus we should be careful in interpreting the observational data 
to compare with theoretical predictions of various cosmological models. 

In this Letter we present a simple analysis of the uncertainty caused by the finiteness of the 
sampled volume, which we call the sample variance. We concentrate on the one-dimensional 
velocity dispersion of the smoothed peculiar velocity field which is the sole quantity to characterize 
its distribution on scales linear theory suffices, provided that primordial fluctuations are Gaussian 
distributed, which we assume throughout the paper. Uncertainties related to the determination 
of the peculiar velocity dispersion include various factors, such as the observational errors in the 
estimation of distances to the objects, sparseness of the number of the clusters/galaxies probed, 
and so on. Monte Carlo calculations using mock catalogues including all these errors have been 
done in the literatures (e.g. Borgani et al. 1997), which are very complicated and require high 
computational costs. In contrast to these approaches, our purpose in this article is to clarify the 
fundamental limitations due to the finiteness of the sample, which is very important in the sense 
that the uncertainty is independent of how accurate we could measure the peculiar velocities of 
galaxies or clusters. 

We develop a simple analytical formula similar to the same kind of analysis on the CMB 
measurement (Scott et al. 1993). We show that the sample variance causes an uncertainty as 
large as 10% on the measurement of velocity dispersion, even if we take a full sphere of radius 
100/i _1 Mpc around us as the observational volume. This is comparable to the current observational 
error (Bahcall &Oh 1996). 



2. Formulation 

First we define the peculiar velocity field, V(x), smoothed over a radius R by 

V R (x) = J V(x)W R (x - x')d 3 x', (1) 

where Wr(x — x') is the Gaussian window function given by W R (x) = (27r)- 3 / 2 ir 3 exp (-x 2 /2R 2 ). 
We define the line-of-sight component by 

V R {x) = Vr(x) ■ (2) 
| a; | 

It is true that, in reality, we must first project and then smooth, but we ignore this incommutability 
since we consider a survey region whose dimension is much larger than R. Otherwise the statistical 
average would become nonsense. To be specific we assume that the observed region is a sphere 
with radius L and denote its volume by Vl = f 71 "-^ 3 with L 3> R. 

One dimensional velocity dispersion estimated from the line-of-sight component in the surveyed 
volume Vl is written as 

X(L,R) = 1- f V R (x) 2 d 3 x. (3) 

The ensemble average of this dispersion is of course independent of the size of the observed volume, 
and it is simply the one-dimensional velocity dispersion with smoothing radius R, which is denoted 
as 

(X(L,R))=a 2 1D (R). (4) 
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Here and hereafter, the angular bracket represents an average over ensembles of the "universes". 
The error associated with the estimation of X(L, R) due to the finiteness of the sampled volume is 
characterized by the variance, 

[X(L, R) - (X(L, R))] 2 ) = (x 2 (L, R)) - (X(L, R)) 2 . (5) 

We assume that primordial density fluctuation is isotropic random Gaussian, and apply linear 
theory, since we adopt large enough smoothing scales suggested by Bahcall et al. (1994). Using 
properties of the Gaussian distribution, we obtain 

(X 2 (L,R)) = ^ J Sx J <Py(v R (x) 2 V R (y) 2 ) 

= ^ J d 3 x J d 3 y ((VrOc) 2 ) (vfc(y) 2 ) + 2 (VH(x)VR(y)) 2 ) . 

Hence the variance of X(L, R) reads 

(X 2 (L,R)) - (X(L,R)) 2 = ^ jf d 3 xj v d 3 y(y R (x)V R (y)) 2 . (6) 

Here (v R (x)V R (y)^ depends only on the geometry decided by x and y, namely, on the three 
variables \x\ = n, \y\ = r 2 and the angle 6 between x and y. Let us normalize the above dispersion 
by (X(L,R)) 2 = af D (R) = (v R (x)V R (x)f to define 



(X(L,R)) Vl Jo Jo Jo 



(X(L,R)) 

with 



V R (x)V R (y) 

Hrur 2 ,e)= . (8) 

y R (x)V R (x)) 



Next we calculate the correlation function f(r\,r 2 ,9). After some algebra, we can relate it with 
parallel and perpendicular velocity covariance functions, $|| and $j_, defined by Gorski (1988) as 

f(ri, r 2 ,9) = $||(ri 2 )(n 2 cos 9 + n x n y sin 6) + ®±(r 12 ) [cos 9n 2 - sin 9n x n y } , (9) 

where, r± 2 , n x and , n y are defined by 

2 2 2 2 n r\cos9 — r 2 rising 

r 12 = \x-y\ =r 1 +r 2 -2rir 2 cos9, n x = , n v = . (10) 

ri2 r 12 

$11 (r) and $j_(r) are the velocity correlation functions normalized to unity at r = 0, or 
$H (0) = $j_(0) = 1, and they are written in terms of the power spectrum of linear density 
fluctuations, P(k), as follows (Gorski 1988). 

$H (r) oc J™dkP(k)W R (k) 2 (j (kr)-2 J -^y (11) 

$ ± (r) oc J™ dkP{k)W R {kf 3 -^±. (12) 
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In the above expressions Wr{U) = exp(— k 2 R 2 /2) is the Fourier transform of Wij(a;), and j m (z) is 

the spherical Bessel functions of the m— th order. 

In the same manner we can also calculate the sample variance of density fluctuations smoothed 

over a radius R at the level of the linear theory. In this case we only have to replace the line-of-sight 

peculiar velocity Vr{x) by the density contrast smoothed over the same radius R, 5r(x). Thus 
VR(x)VR(y)\ in equation is replaced by (5R(x)dR(y)), and / in equation @ is replaced by 
'(^12) = C_R( r i2)/C/?(0), where £r(t) is the two point correlation function of the smoothed density 

field. 

In this article, we adopt the power spectrum of cold-dark-matter (CDM) models given in 
Efstathiou, Bond & White (1992) as 

-2/7* 



P(k) = Bk{l+\uk + {f3kfl 2 + ( 1 k) 2 Y) * (13) 



where a = (6.4/r)/ i - 1 Mpc, (3 = (3.0/r)/i^ 1 Mpc,7 = (U/Tj/T^Mpc, p = 1.13, and B is the 
normalization factor. Since we are dealing with normalized quantities in the linear theory, the 
system is characterized by the shape parameter V = QqJi, and the smoothing radius R only. 



3. Results 

First we plot $ii(r),$j_(r), and H(r) in Figure 1. Apparently, the peculiar velocity field has 
larger correlation length than the density field, in particular, $j_(r) has a very large correlation 
length. This is understandable because in linear theory the density contrast, 5, and the peculiar 
velocity, V, satisfy the relation V cx V • V~ 2 <5, which implies that the power spectrum of the 
peculiar velocity field, P v (k), is related to P(k) as P v (k) oc k~ 2 P(k) (Peebles 1980). Therefore, 
velocity field is more weighted to smaller k and hence has a much larger correlation length, which 
causes a larger sample variance as we see below. 

In Figure 2, we plot the sample variance of the one-dimensional velocity dispersion, £' 1 / 2 (L, R), 
with smoothing radius R = 3(0.2/r)/i" 1 Mpc and 6(0.2/r)/i _1 Mpc. For comparison we also plot the 
corresponding sample variance for the density field in linear theory. To simplify our presentation, 
the length scale in the horizontal axis is shown in unit of (0.2/r)/i _1 Mpc. 

Apparently, the error due to the finiteness of the observed volume is much larger for the 
velocity field than for the linear density field Q. This figure shows clearly that even with a spherical 
observed region with radius 100/i -1 Mpc, the error on the estimation of the one-dimensional velocity 
dispersion, which is approximately equal to E l l 2 (L,R)/2, is as large as 10%. This is comparable to 
the observational error of ~ 10% in Bahcall and Oh (1996). Note that magnitude of the error we 
have obtained here is for an ideal case that a specific region was observed fully and homogeneously, 
and that the actual sample variance would be even larger with the same depth of the sample. 



Note that linear theory is insufficient to reproduce the magnitude of density fluctuations on the scales we are dealing with. 
Hence we do not claim the above result reflects the correct sample variance of the density field. All we would like to stress here is 
that magnitude of the sample variance is sensitive to the functional shape of the power spectrum. 
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4. Effects of the peculiar velocity of the observer 

So far, we have not specified the peculiar velocity of the observer, who is supposedly at the 
center of the surveyed sphere with radius L. The peculiar velocity of the observer can be estimated 
from the dipole anisotropy of the cosmic microwave background radiation. If we assume that 
it is totally due to our peculiar motion then COBE observation gives V Q b s = 627km/s (Kogut 
et al. 1993). One may suspect that if the correlation of the peculiar velocity field is so strong, 
the observed velocity dispersion in the sphere may strongly be affected by the value of V b s . 
Here we estimate how much the observed velocity dispersion is expected to shift as a function of 
v = V obs /a 1D (R). 

The conditional probability distribution function of the radial peculiar velocity has been 
calculated in Gorski (1988, eq. [7]) with a specific value of the peculiar velocity of the observer, 
from which we find 

[V£(r)) =a 2 1D (R)\l-^(rf(l-u 2 co S 2 e)], (14) 



where (■■■)„ denotes the ensemble average on condition that the observer has a peculiar velocity of 
\Vobs\ = v(Tid(R), and is the angle between V \, s and r. We thus find the constrained ensemble 
average of X(L,R) is given by 



(X(L,R)) u = ^l f dMl-$||W 2 (l-- 2 cos 2 6)}. (15) 



<tid{R)\ 

Therefore the fractional change of the average due to the additional constraint is given by 

(X(L, R)) v - (X(L, R)) v 2 -3 f L ^ , , 2 _ 



(X(L, R)} I? 



f dr$ii(r) 2 r 2 . (16) 
Jo 



For COBE-normalized CDM model (with h = 0.75, T = 0.2) we find v = 1.58 and the above 
change is only about 0.5% for R = 3(0.2/r)/ i - 1 Mpc and L = 100(0. 2/r)/i- 1 Mpc, much smaller 
than the sample variance we have discussed. Hence the peculiar motion of the observer does not 
cause a serious problem. 



5. Discussion 

We have shown that the fact that peculiar velocity field is very sensitive to small wave-number 
modes implies that we must observe a sufficiently large region in order to extract useful information 
on cosmological models from observational data. We have dealt with a smoothed peculiar velocity 
field, taking the smoothing length large enough to warrant the validity of linear theory in which 
the only important quantity is the dispersion as long as the primordial fluctuations are distributed 
Gaussian. Hence we concentrated on the sample variance on the dispersion of the one-dimensional 
line-of-sight peculiar velocity and have shown that it is not negligible at all. 

Although estimation of cosmological parameters from the analysis of peculiar velocity field 
mostly uses more contrived methods {e.g. Dekel 1994), it would be natural to expect that similar 
limitation applies in these approaches as well, because they essentially start with the same kind of 
observational measures, namely, the line-of-sight peculiar velocity, and because linear theory has 
been shown to suffice on the relevant scales (Bahcall et al. 1994). In order to give a rough idea of 
magnitude of the error on the estimation of cosmological parameters caused by the sample variance, 
let us take the smoothing radius large enough that the linear theory applies for the density contrast 
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as well, say R = 12h Mpc, and relate the root-mean-square density fluctuation obtained from 
galaxy distribution, a ga i(R), with the one-dimensional peculiar velocity dispersion <jid(R). Using 
the formula of Colberg et al. (1997), we find a\ D {R) = %- 2 a ga i{R) 2 F (V ,R)/b 2 , where b is the 
bias parameter and F(T, R) is a function decided only by the shape of the matter power spectrum 
and R. Even in the ideal and hypothetical case that both a ga i{R) and the coefficient F(T,R) are 
known exactly, the above simple formula tells us that we have an inevitable uncertainty in 0°' 6 /6 
of about 25% for L = 60/i _1 Mpc and 20% for L = lOO/i^Mpc (in the case CDM spectrum flT|) 
with T = 0.2 is adopted). 

In conclusion, the strong dependence of the peculiar velocity field on small wave-number modes 
has the two aspects; the advantage is that the linear theory suffices even on relatively small scales, 
while the disadvantage is that we must have a larger observational volume to measure it accurately. 
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Fig. 1. — The normalized correlation functions of velocity and density with smoothing radius 3(0.2/r)/i _1 Mpc. 
The solid curve represents the two point correlation function of the linear density field, £(r). The dotted 
line represents the parallel velocity covariance function, $ii(j"), and dashed line represents the perpendicular 
component, 3>j_(r). 




Fig. 2. — E l l 2 (L,R), the expected error on the estimation of the universal velocity dispersion. The horizontal 
axis is the radius L of the observed sphere in unit of (0.2/r)/i _1 Mpc. The dotted line corresponds to E l l 2 (L, 6), 
and the solid line to -E 1//2 (L,3), both for the one dimensional peculiar velocity. 

The long-dashed and short-dashed lines represent the corresponding quantities for the linear density contrast. 



